Passive Scalar Structures in Supersonic Turbulence 
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We conduct a systematic numerical study of passive scalar structures in supersonic turbulent flows. We find 
that the degree of intermittency in the scalar structures increases only slightly as the flow changes from transonic 
to highly supersonic, while the velocity structures become significantly more intermittent. This difference is due 
to the absence of shock-like discontinuities in the scalar field. The structure functions of the scalar field are well 
described by the intermittency model of She and Leveque [Phys. Rev. Lett. 72, 336 (1994)], and the most intense 
scalar structures are found to be sheet-like at all Mach numbers. 
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PACS numbers: 

Supersonic turbulent motions have been observed over a 
wide range of length scales in the interstellar medium JT1, 
and the mixing of heavy elements released from stellar winds 
and supernovae occurs in such a supersonic turbulent envi- 
ronment. Thus understanding passive scalar physics in super- 
sonic turbulence is crucial for the interpretation of many ob- 
servational results concerning cosmic chemical abundances. 
The scalar field we study here is the concentration field, 
0(x, t), of passive tracers. 

The statistical approach commonly used to study struc- 
tures in turbulent systems is to analyze structure functions 
(SFs), defined as Sp(r) = (\8v(r)\ p ) for the velocity field, or 
5|(r) = (\56(r)\P) for the scalar field. Here Sv(r) and 86{r) 
are, respectively, the longitudinal velocity increment and the 
scalar increment over a distance r. In the scale range where 
the non-linear interactions prevail, the SFs are expected to fol- 
low power laws, Sp(r) oc r^ v ^ p ', and Sp(r) oc r^ 8 ^. 

The velocity scaling exponents £ v (p) m incompressible tur- 
bulence depart significantly from the prediction, £ v (p) =p/3, 
by the similarity hypothesis of Kolmogorov at high orders 
(p > 3), a phenomenon known as anomalous scaling or in- 
termittency [2]. The physical origin of the departure is the 
strong spatial fluctuations of the energy dissipation rate, which 
provides intermittency corrections to £ v (p)- Perhaps the most 
successful intermittency model is the one by She and Leveque 
(hereafter the SL model) [3], and its prediction for £ v (p) m 
incompressible turbulence agrees with experimental results at 
high accuracy [4]. 

In supersonic turbulence, the existence of velocity shocks 
gives more intermittent velocity scalings, and the degree of 
intermittency increases with the Mach number, M. Padoan et 
al. showed that the velocity scaling in supersonic flows can 
be unified using the SL model with a parameter, the dimension 
of the most intense velocity structures (MIVSs), that varies 
with M. These structures make a transition from filamentary 
(vortex tubes) at small M to sheet-like (shocks) at large M. 

Studies of passive scalars in incompressible turbulence 
found that scalar SFs also have anomalous scaling behav- 
iors, and that scalar structures are more intermittent than the 
velocity structures |@]. The most intense scalar structures 
(MISSs) are sheet-like, corresponding to the observed cliffs 
in the scalar field iM- llOll . In this Rapid Communication, we 
conduct a systematic numerical study of passive scalar struc- 



tures in supersonic turbulence. 



The classic theory for turbulent mixing Hi 111 assumes 
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where eg is the average scalar dissipation rate. Refining eq. 
((TJ to account for the fluctuations in the scalar dissipation rate 
gives 89(r) 2 ~ e$(r) S J^ where eg(r) is the dissipation rate 
averaged over a scale r. The quantity, 8v(r)S9(r) 2 , associ- 
ated with the scalar cascade is of special interest [9']. Defin- 
ing "mixed" SFs as S£(r) = (\5v(r)S8(r) 2 \ p / 3 ), the inter- 
mittency corrections to the mixed structures come completely 
from the fluctuations in the scalar dissipation rate. We study 
the velocity, the mixed, and the scalar SFs, denoted as Sp(r), 
S^(r) and Sp(r), respectively. We use super- or subscripts, 
v, m and 6 to specify each case, and when no such super- or 
sub- scripts are used, the discussion is general, referring to all 
the three cases. 

We use the 512 3 simulation data from Pan and Scanna- 
pieco lll2"ll . who simulated isothermal hydrodynamic turbulent 
flows at 6 Mach numbers (M = 0.9, 1.4^2.1, 3.0, 4.6, and 6.1) 
using the FLASH code (version 3.2) lfl3ll . We integrated the 
advection equation for the concentration field, 



ft. 



-v ■ V6> = S(x,t), 



(2) 



where S(x, t) represents the tracer sources. In each flow we 
evolved three independent scalars with statistically equivalent 
source terms. Averaging over the three scalars gives more ac- 
curate measurements. Theoretical studies predicted the exis- 
tence of an inverse scalar cascade in highly compressible tur- 
bulence 11411 . which, however, was not observed in our simu- 
lations, as the critical compressibility for the inverse cascade 
is not reached even at M = 6.1. Several lines of evidence for 
a direct scalar cascade at all M are given in Ref. [12]. 

The viscous and diffusion terms are not explicitly included 
in our simulations, and the dissipation of kinetic energy and 
the scalar variance is through numerical diffusion. With the 
dissipation rates determined from their balance with the forc- 
ing rates, we estimated the effective viscosity and diffusivity 
using the relations between the dissipation rates and the veloc- 
ity/scalar gradients. The Taylor Reynolds number is estimated 
to be ~ 250 in the M = 0.9 flow. It decreases to ~ 140 at 
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M = 6.1 because the code induces a larger effective viscosity 
to stabilize stronger shocks. One important issue in passive 
scalar physics is the effect of the Schmidt number, Sc (see 
Ref. 11511 for a detailed review). Our estimated viscosity and 
diffusivity give an effective Sc ~ 1 at all M. However, we 
point out that uncertainty exists in the estimate, and the Sc 
effect in our simulations cannot be exactly evaluated. The Sc 
effect for mixing in supersonic turbulence needs to be studied 
in future simulations including viscous and diffusion terms. 

We first consider the 3rd order SFs, which will be taken as 
references in our scaling analysis below. The measured values 
for £ v (3), £m(3) and &(3) are listed in Table 1. At M = 0.9, 
we find that £ v (3) and £ m (3) are close to unity, consistent with 
Kolmogorov's 4/5 law [2], and the 4/3 law by Yaglom fl^l in 
incompressible turbulence. As M increases to 2, £ v (3) for the 
velocity field increases, and then saturates at 1.22 for M 3, 
close to that found in Ref. lfl7ll . 

Eq. (Q} suggests that £ m (3) is equal to 1 at all Mach num- 
bers. In fact, we find that £ m (3) is near 1 at each M, support- 
ing the general validity of the classic theory for passive scalar 
physics in supersonic turbulence. However, we note that £ m 
starts to increase continuously at M J; 3, meaning that eq. ([TJ 
is not exactly obeyed at large M. This is likely to be caused by 
the effect of compressible modes on scalar structures, which 
is not reflected in eq. ([T). 

For scalar structures, we find £e(3) = 0.87 at M = 0.9, 
which is the same as that obtained in Ref. for the incom- 
pressible case. The exponent £g(3) first decreases with in- 
creasing M, and then starts to increase at M > 3, similar to 
the behavior of £g(2) found in Ref. 1 12]. A detailed explana- 
tion for this behavior is given in Ref. Ill 21 . 

Measuring scaling exponents £(p) at high orders is difficult 
due to the limited inertial range in our simulations. Here we 
exploit the concept of extended self-similarity 11811 . Plotting 
the SFs at all orders against the 3rd order ones gives more 
extended power-law ranges, allowing more accurate measure- 
ments. We thus measure the scaling exponents, relative 



to the 3rd order SFs, defined as S p (r) oc (S^r)) 



C(p) 



By 



definition, £(p) = f(p)/f(3) and C(3) = 1. 

Fig. 1 shows the measurement of £(p) for the scalar SFs 
at M — 6. The data points in the figure, from left to right, 
correspond to r = 4 to 196 (in units of the size of the com- 
putation cells). An extended power-law range is seen at low 
to intermediate orders, and the range extends into the dissipa- 
tive scales [18]. At large p, the power-law range becomes 
smaller, and we calculated the slope using least-square fits 
with the central 6 points, corresponding to r 6 [12, 64]. An 
examination of the probability distribution of Sv(r) and S9(r) 
shows that our data have good statistics at orders p up to 10 
for M & 3. At M = 4.5, and 6.1, the statistics forp ~ 10 
appears to be sufficient, but the measurement uncertainty be- 
comes larger, especially for the velocity SFs. One reason for 
this is that, at larger M, the effective numerical viscosity in 
our simulations is larger, giving a narrower inertial range . 

In Fig. 2, we plot the measured exponents, £(p), for the 
three kinds of SFs at M = 0.9 and M = 6.1. In the M = 0.9 
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FIG. 1: Compensated scalar SFs in the M = 6.1 flow at orders, p, 
from 1 (bottom) to 10 (top). Horizontal lines show the quality of the 
power-law fits. For clarity, SFs at p — 1 and 2 are multiplied by 0.3 
and 0.7, respectively. 



flow, ((p) for all the three SFs agrees well with the results 
for incompressible turbulence (within 4% at all orders, p) 
JSSEllli. At M = 0.9, ( 6 {p) is smaller than both ( m (p) 
and Cv(p)i indicating the scalar structures are the most inter- 
mittent in the incompressible limit. As M increases to 6.1, 
the relative degree of intermittency is reversed, with the ve- 
locity field becoming more intermittent than the scalar field. 
The transition between these two regimes occurs at M = 2, 
where the £ v (p) and £e (p) curves are very close to each other. 

The scaling exponents for the velocity SFs decrease steadily 
and significantly with increasing M. This is because the fre- 
quency and the strength of shocks in supersonic flows increase 
with M. Shocks are strong intermittent structures, which tend 
to decrease the scaling exponents at high orders, as illustrated 
by Burgers turbulence • 

In contrast, there is no such effect for the scaling expo- 
nents of the scalar structures, which decrease only slightly 
as M changes from 0.9 to 6.1. No discontinuous structures 
like the velocity shocks are expected in the scalar field at any 
AI. As discussed in Ref. 11211 . velocity shocks can squeeze 
the scalar field, and amplify the concentration gradients, lead- 
ing to edge-like structures around the shocks (see Fig. 1 in 
Ref. 112J). However, these edges are not discontinuities be- 
cause the amplification factor for the scalar gradient by a 
shock is finite, about equal to the density jump across the 
shock. The continuity in the scalar field is also expected 
from the conservation of tracer mass. The tracer density jump 
across a shock is the same as that of the flow density, and thus 
the concentration, the ratio of the two densities, is continuous 
across the shock. The lack of shock-like discontinuities in the 
scalar field is responsible for the different scaling behavior of 
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FIG. 2: Scaling exponents £(p) for SFs at A/ = 0.9 and M = 6.1. 
Note that <(p) = £(p)/f (3) and C(p) / f (p) when £(3) / 1. Lines 
are predictions of the SL model, and parameters used in the model 
are given in Table 1 . 



TABLE I: Scaling exponent £(3) of 3rd order structure SFs, inter- 
mittency parameter /3 and the fractal dimension d of MISs 
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d v 
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0.96 
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the scalar SFs as a function of Al from the velocity SFs. 

We next show that the scaling behaviors of the scalar struc- 
tures in supersonic turbulence can be well fit by the SL in- 
termittency model yfl, which has achieved substantial success 
in a wide range of turbulent systems 10, llol] . The model is 
based on the consideration of a hierarchy of structures, charac- 
terized by ratios of successive SFs, F p (r) = S p +i(r) / S p (r). 
With increasing p, F p (r) represents structures of higher inten- 
sity level, and (r) corresponds to the most intense struc- 
tures (MISs) at the scale r. The main assumption of the model 
is the existence of a symmetry in the hierarchy, 



F p+1 (r)=A p F p (rfF 00 (r) 1 -P, 



(3) 



where A p is assumed to be independent of r, and (3 is referred 
to as the intermittency parameter. 

Solving the recursion relation, eq. ((3), gives, 



FIG. 3: Testing symmetry eq. ((3} in the M = 6.1 flow. For clarity, 
data points for mixed and scalar structures are shifted upward by 0.3 
and 0.6, respectively. 



where 7 is the scaling exponent for the MISs, i.e., F^fr) oc 
r 7 . The parameter C has a physical interpretation as the codi- 
mension of the MISs, defined as C — 3 — d in 3D with d 
being the fractal dimension of the MISs. An appealing feature 
of the SL model is that, through the assumed symmetry, the 
entire hierarchy of structures is related to the MISs. Therefore, 
if applicable, the model provides an approach to measure the 
physical properties (7 and C) of the MISs. 

If eq. (0]l applies for £,(p), a similar equation exists for ((p), 



c(p)=7 / P+c(i-/n 



(5) 



where 7' = 7 /£(3) and C = C/f(3). Since C(3) = 1 by 
definition, we have C = (I - - /3 3 ). We will deter- 

mine (3 and C from our simulation data. Unlike 7 and C, the 
physical meanings of 7' and C are not straightforward when 
£ (3) ^ 1. Therefore, after obtaining 7' and C", we convert 
them to 7 and C using the measured values of £(3). In par- 
ticular, we will calculate dby d = 3 — C = 3 — £(3)C". We 
note that the dimension derived in Ref. [5] for MIVSs in su- 
personic flows was defined as 3 — C. It is different from our 
definition of d and will be referred to as d! . 



Eq. (0) can be rewritten as 112 ill . 



F P+ i(r) 
F 2 (r) 



Ml 

Fi(r) 



(6) 



We can thus examine the validity of the assumed symme- 
try, and hence the applicability of the SL model, by check- 
ing whether the \og(F p+ i / F2) - \og(F p /Fi) curves have the 



£(p) = 7P +c(i-/n 



(4) 



same slope at all orders, p 12111 . In particular, if A p is inde- 
pendent of p, these curves would fall onto the same line. This 
is indeed the case for the velocity structures in incompressible 
flows [22]. The same is found here for the velocity structures 
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at M = 0.9, 1.4 and 2. In these three cases, we estimated 
(3 V by fitting a single line to the data points from all orders, 

In all other cases, A p is not constant, and the log(F p+1 /F 2 ) 
- log(F p /Fi) curves at different p are not on the same line. 
However, these curves are generally parallel to each other, 
also confirming the validity of eq. ([3J. The symmetry test 
for M = 6.1 is shown in Fig. 3. Similar results are found 
for other M. The 6 data points at each order correspond to 
the scale range r € [12, 64]. Small variations are observed in 
the slopes at different p, and we evaluated j3 by averaging the 
slopes over orders up to p = 10. 

The symmetry test determines the parameter j3, and with 
0, we obtain C by fitting eq. Q to the measured exponents, 
C(p). The dimension, d, is then calculated by d = 3 — £(3)C". 
The results for M = 0.9 and 6.1 are shown in Fig. 2. A good 
match is seen between the model predictions (lines) and nu- 
merical results (data points). The measured values of f3 and 
d are listed in Table 1 . Note that, with increasing M , /3 V de- 
creases significantly, while the change in fig is slight. This 
corresponds to the trend of the scaling exponents observed 
in Fig. 2. The parameters obtained in our M = 0.9 flow 
agree well with those from studies for incompressible turbu- 
lence fl SHHl. Consistent with Ref. H, we find that, with 
increasing M, the MIVSs make the transition from filamen- 
tary to sheet-like, with d v changing from 1 to 1.8 as M goes 
from 0.9 to 6. 1 . This range of d v corresponds to the change of 
d' v from 1 to 2, in agreement with d' v measured in Ref. ||5|]. 

Of particular interest here is the dimension, dg, for scalars. 
Extensive experimental and numerical results have shown that 
the MISSs in incompressible flows are in the form of sheets, 
known as the cliff structureslM, 12311. and the dimension dg ~ 2 
derived in our M = 0.9 flow is consistent with these results. 
Furthermore, the MISSs have a dimension of ~ 2 in all our su- 
personic flows. This result can be understood by considering 
the possible effects of compressible modes in shaping the ge- 
ometry of the intense scalar structures. As mentioned earlier, 
strong compressions in supersonic flows can produce edge- 
like scalar structures across shocks [12]. The strong shear that 
usually exists behind shocks can further make scalar sheets in 
the narrow postshock regions. These may provide the main 
contributions to the MISSs. Although formed from a different 
mechanism, they are geometrically similar to the cliff struc- 
tures found in incompressible flows. Therefore whether or not 
their contribution to the MISSs is dominant, dg is expected to 
be ~ 2 at all M. 

The dimension, d m , for the mixed structures shows a more 
complicated dependence on M than the velocity and scalar 
structures. As M increases from 0.9 and 2, it decreases from 
1.9 to 1.5, while d v increases and dg is essentially constant. At 
larger M, d m increases again, and the MISs for all the three 
cases have d close to 2 . 

In conclusion, we carried out a systematic study for passive 
scalars in supersonic turbulence. We find that, with increasing 



M, the velocity scalings becomes significantly more intermit- 
tent, but the degree of intermittency in the scalar structures 
increases only slightly with M, This is because, at any Mach 
number, the scalar field does not have discontinuous structures 
like velocity shocks. The scalar scalings at all Mach num- 
bers are well described by the intermittency model of She and 
Leveque. Fitting the model prediction to the measured scal- 
ing exponents shows that the most intense scalar structures are 
sheet-like at all Mach numbers. 
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